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Abstract 

A systematic study of the effects of polarization mode dispersion on broad-band and narrow-band, single pump, 
fiber parametric amplifiers is realized through numerical solutions of the equations governing the interaction. The 
nonlinear polarization rotation is shown to be a relevant effect that can increase gain randomness when it mixes with 
polarization mode dispersion. In unidirectionally spun fibers the signal-pump alignment can be highly increased and 
the gain enhanced. However, in spite of the enhanced alignment, large polarization mode dispersion phase-mismatches 
the interaction and the gain decreases to zero. 

I. Introduction 

Fiber optical parametric amplifiers (FOPAs) have been under investigation for many years. In the single pump 
configuration both broad-band (BB-FOPA) [1] and narrow -band (NB-FOPA) [2], [3] amplification can be obtained. 
The former has been mainly addressed for amplification and wavelength conversion [4] while the latter yields a 
superb technique for slow and fast light generation [5], [6]. 

FOPAs are highly affected by the fiber random birefringence (polarization mode dispersion - PMD) [7], [8]. In 
an ideal, perfectly isotropic fiber the waves would maintain their input states of polarization (SOP) over the entire 
fiber, the gain would be deterministic and it would depend only on the initial SOPs [9]. In low birefringence fibers 
the evolution of the pump, signal and idler SOPs is random, due to the PMD. Eventually, the gain is reduced and 
becomes a random quantity, as theoretically analyzed in [7] for BB-FOPA. 

Recently, the benefits of unidirectional fiber spinning in Raman [10] and Brillouin fiber amplifiers [11] have been 
underlined. Here, the same detailed study is extended to single pump, BB- and NB-FOPAs. Note that, differently 
from the Raman and Brillouin fiber amplifiers, in FOPAs pump polarization scrambling cannot be used to avoid 
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polarization dependence. Other techniques like dual pumping [12], [13] or polarization diversity [14], [15] were 
proposed to overcome the problem, but at the cost of an increased set-up complexity. So, the results presented here 
are particularly relevant because unidirectionally spun fibers can represent a breakthrough for implementing highly 
performant BB- and NB-FOPA with simple experimental setups. 

II. Propagation model 

The parametric interaction of the pump, signal and idler monochromatic complex Jones vectors, \Ap), \ As), \ Ai), 
in the undepleted pump approximation and neglecting the nonlinear effects of the signal and idler on the pump, is 
governed by [7]: 

dz\Ap) = j {jap + k{ojp) - 0{ojp) ■ a/2) + 
+.n/3 {2{Ap\Ap) + \A;){A;\) \Ap), 
dz\Ass) j {jas,i + k{us,i) ~ P{^s,i) ■ cfji) + 
+J27/3 {2{Ap\Ap) + \Ap){Ap\ + |^;)(^;|) lA,,,) + 

+J7/3 {{Al\Ap) + 2\Ap){Al^ (1) 

where 7 = 2 W^^kmT^ is the nonlinear coefficient, , h~ p, s, i are the loss coefficients at the optical pulsations 
of the pump, signal and idler LOp,ijJs,uji that satisfy simultaneously the energy 2ijjp = ujs + oji and the nonlinear 
phase matching condition Afc = fc(ijJs) + k{u!i) — 2k{u!p) w k2p {uJs — '-^pY + ^ip {'^s — ^^p)*/12 = ~2^Pq, {knh is 
the n-th angular frequency derivative of k{u!), evaluated at uj^ [1] and Pq is the pump power). Transition from BB- 
to NB-FOPA can be obtained by tuning ujp in the vicinity of the fiber zero dispersion wavelength, thus modifying the 
dispersion coefficient k2p. When the sign of k2p is negative, gain occurs at uJs within a broad band around the pump 
frequency. For positive k2p and if k4p is negative, the nonlinear matching can be satisfied in two narrow spectral 
bands, largely detuned from the pump pulsation. The following parameters have been used in the simulations: 
BB-FOPA, k2p = -1.18 X lO^^s sVm, k^p = 1 x IQ-^^ s^/m [7];NB-FOPA, fcsp = 7.68 x lO^^s sVm, 
kip = -5 X 10-55 s4/jj^ [3]^ 

The Stokes vector f3{z,uj) is the local birefringence vector and a is the vector of the Pauli spin matrices [19]. 
For unspun fibers, (3{z,uj) = can be obtained through the so called random modulus model (RMM) 

[16], that describes the PMD phenomenon and effects through three characteristic lengths: the beat length (Lb), the 
birefringence correlation length (Lp) and the polarization correlation length (Lc)- The first depends on the pulsation 
through Lsii^) = cjo^bI'^o)/'^; the latter is independent of lo {Lp = 9 m was fixed in the numerical integrations). 
The first two contribute to determine the PMD coefficient (PMDC), hereinafter defined as D = {/S.t'^) / L where 
(Ar^) is the fiber PMD mean square differential group delay and L(= 2 km) is the fiber length. The polarization 
correlation length Lc, is the length scale over which the SOP is modified and depends on the previously defined 
lengths and on the fiber spinning properties. This length scale is a fundamental parameter when nonlinear PMD 
sets in [16], as in this case. 
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For unidirectionally spun fibers the birefringence vector becomes (3{z,uj) ~ i?3[20(z)]/3„s(z, w) where is a 
Mueller matrix representing a rotation around the third component of the Stokes space (U3). The angle of rotation is 
given by the spin function <f>{z) = 2ttz/p where p is called the spin pitch. It was found that unidirectional spinning 
[20]: a) decreases the PMDC and b) modifies the polarization correlation length Lc, according to: 



aiT^Lp \ p' 



for L\» Att'^L'^,/ + (167r2L2^)/p2 and where a = 4 (a = 2) for linear (circular) input SOPs. 

A few thousands statistical realizations of the stochastic processes, using the RMM, and subsequent integrations 
of Eqs. [T] have been performed to determine the mean gain /iq = {{As{L)\A*{L))) / {As{0)\A*{Q)) and the 



relative gain standard deviation (STD) ctg ^ y {{As{L)\Al{L))) / {^^ {As{L)\Al{L)))'^ - 1. Numerical results are 
represented in all figures as functions of the PMDC that is calculated from the statistical ensemble. 

III. Numerical results 

For unspun fibers, the mean gain and STD from numerical solutions are presented (circles) for initially linear and 
parallel SOPs in Fig. [Ufor BB-FOPAs (A^ ~ 1550.2 nm. A,. ~ 1514 nm) and in Fig. |2]for NB-FOPA (A^ ~ 1541 nm. 
As c± 1387 nm) as a function of D. Mean and STD deviation for other input SOPs are not presented here, for the 
sake of brevity, but all the results show the same behavior that is described below. When the PMDC is large, 
the gain decreases down to zero, because PMD tends to phase mismatch the waves [7], [8]. In this regime the 
predictions of the mean Eqs. of [7] (solid curves in both figures), where nonlinearity was averaged and an effective 
nonlinear coefficient 7e = 87/9 was found [17], [18], are recovered both for BB- and NB-FOPA. The PMD affects 
more NB-FOPA than BB-FOPA mean gain; this can be explained by the fact that the effects of PMD, at the same 
value of D, increases with the frequency detuning among the waves. 

However, for small PMDC the mean and STD values do not correspond to the theoretical findings. The mean 
gain is actually larger and the STD presents a peak (see ctg in Figs. [T] and |2] for D « 10^'^ ps/Vkm) not predicted 
by the theory. The discrepancy stems from the nonlinear PMD [16], induced by the nonlinear polarization rotation 
(NPR) effect, as already briefly described in the case of Raman amplification [21]. In the Stokes space, the pump 
vector P = {Ap\(j\Ap), is governed by: 

d,P = [/3(z, Up) + 27/3(P • U3)U3] X P. (3) 

When the pump polarization evolution length Lc is close to the NPR length scale NPR becomes a significant 
source of randomness [16], which finally leads to enhance signal gain randomness. To prove this affirmation, the 
theoretical values of Lb satisfying the condition Lc = Ln pr in Eq. |2l in function of the spin rate v ~ 1/p, have 
been calculated for the BB-FOPA and are shown in Fig. [3] (solid curves), in the limit of validity of Eq. |2] (similar 
results are obtained for the NB-FOPA). The value of the NPR length was estimated as Lnpr ~ 9/{2^Po) = 2.25 
km in the hypothesis, which is valid if the nonlinear mixing is strong enough, that the mean power in the stochastic 
vector (P • ^3)^3 is one third of the total pump power. In the same Fig. [3] the values of for which the STD 
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reaches the relative maximum in the numerical solutions are also shown for comparison (circles). The agreement 
is very good both for linear (black solid curve and circles) and circular (grey solid curve and circles) input SOPs. 

NPR effects become relevant for rapidly spun fibers, for which the averaging of nonlinearity cannot be performed 
and the effective nonlinear coefficient is 7 [22]. In fact, as the spin rate increases, the peak in the STD, due to the 
NPR, shifts to larger PMDC values because D does not follow the same scaling law with p of Lc [20]. Moreover, 
the mean gain of both BB- and NB-FOPA increases because of the enhancement in the mean pump-signal SOPs 
alignment {p- s) = {cos{9)), as already demonstrated in [10], [11]. This effect is shown in Fig. |4]for the NB-FOPA. 
Finally, observe that highly increasing the spin rate does not lead to a gain enhancement similar to that observed for 
Raman and Brillouin amplifiers [10], [11]. In fact, PMD affects the pump-signal alignment, restored by spinning, 
but also the phase matching [7], [8], whose dependence on spinning is actually neglectable. 

IV. Conclusions 

In conclusion, the effects of polarization mode dispersion and unidirectional fiber spinning in broad-band and 
narrow-band single pump, fiber optical parametric amplifiers were studied by means of numerical solutions of the 
governing equations. The effects of nonlinear polarization rotation are relevant in such amplifiers, in particular when 
unidirectional spinning is applied. The gain can be enhanced by spinning and, in NB-FOPA for 10~^ ps/\/km < 
D < 4 ■ 10~^ ps/Vkm, the gain variance can be decreased. This observation is relevant for slow and fast light 
generation, where polarization mode dispersion highly affects the interaction [6], [8]. The gain increase stems from 
the spinning enhancement of the alignment between the pump and the signal along the fiber However, differently 
from Raman and Brillouin amplifiers, the phase matching condition is also affected by polarization mode dispersion 
and the fiber spinning is ineffective to counteract mismatch. So, a limit to gain enhancement is found. 
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D[ps X kni-i/2] 

Fig. 1. Mean gain and STD for initially linear and parallel SOPs for BB-FOPA as a function of D. For unspun fibers, the solid curve is the 
result of the mean equations of [7], while the circles are the result of the numerical integration of Eqs.[T] The numerical results for spun fibers 
with spin rate i/ = 1/p = 0.5, 1, 2 and 10 turns/m are represented respectively by diamonds, triangles squares and stars. 
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2: [km] 

Fig. 4. Evolution of (cos{9)) vs. z for in NB-FOPA (D ~ 0.02 ps/\/km). Continuous curves, from top to bottom, refer to spun fibers with 
decreasing spin rate (10, 2, 0.5 tum/m); the dashed curve refer to unspun fibers. The corresponding mean gain are 26.5, 24.9, 22.9, 16.1 dB, 
respectively. 
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